SOME NOTES CONCERNING THE DYNAMICS OF 
NONCOMMUTATIVE LUMPS CORRESPONDING TO 
NONTRIVIAL VACUA IN THE NONCOMMUTATIVE 
YANG MILLS MODELS WHICH ARE PERTURBATIVE 
BRANCHES OF M(ATRIX) THEORY 

o 

O . CORNELIU SOCHICHIU 

Abstract. We consider a pair of noncommutative lumps in the noncommuta- 
, tive Yang— Mills/M(atrix) model. In the case when the lumps are separated by 

■ a finite distance their "polarisations" do not belong to orthogonal subspaces 
' of the Hilbert space. In this case the interaction between lumps is nontrivial. 

We analyse the dynamics arisen due to this interaction in both naive approach 
of rigid lumps and exactly as described by the underlying gauge model. It 

■ appears that the exact description is given in terms of finite matrix mod- 
^ ' els/multidimensional mechanics whose dimensionality depends on the initial 

' conditions. 

I> 

o 
o 

1. Introduction 

Recent progress in theories over noncommutative spaces (for a review see e.g. 
and references therein), is stimulated by their importance for the nonpertur- 
Q^l bative dynamics of string theory 

^ ■ The noncommutative models share some common features with their commuta- 

tive counterparts, however, there is a striking dissimilarity between them in some 
other aspects. One particular feature of noncommutative field theories discovered 
recently and which attracted a considerable interest is that in noncommutative 
models there exists a kind of localised solutions nonexistent in the models on com- 
I mutative spaces. Although they are different from what is a soliton in usual sense, 

these solutions are conventionally called "noncommutative solitons" . In actual work 
we consider a subclass of such configurations. As it appears that the "noncommu- 
tative solitons" in actual work even do not carry energy (at rest) more adequate 
would be the term of "vacuum" or "lump" solution, throughout this paper we will 
keep the last name for them. 

Noncommutative lumps, in a scalar model with a potential having nontrivial 
local minima were first discovered in in the limit of strong noncommutativity. 
These solutions were interpreted as condensed lower dimensional branes living on 



X 



a noncommutative brane ||10|, 11 . They were further generalised to the case of 
a mild noncommutativity by allowing the presence of the nontrivial gauge field 
backgrounds [^-||l^. This solutions correspond to nontrivial gauge field fluxes 
p6[|-[p^. The particular property of lump solutions we are considering in the 
actual paper is that they are "made" purely of the gauge fields. However using the 
equivalence between different noncommutative Yang-Mills-Higgs models [|l9|, pO) , 
these configurations can be mapped into noncommutative solitons in the sense of 
Ref. H or others. 
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The general multi-lump solutions look like sums of projectors to mutually or- 
thogonal finite-dimensional subspaces of the Hilbert space. If subspaces are not 
orthogonal the configuration fails to be a static solution and lumps start to inter- 
act. 

An approach to describe interacting lumps were proposed in ^ by using 
the substitution of the configuration of shifted lumps by a close one but being a 
static solution. In this approach the interaction of the lumps is described by the 
motion in the curved moduli space of static solutions. This approach, however, 
would be valid only provided that the motion is confined to the moduli space of the 
static solution which requires it to be stable. There are, however, indications that 
the noncommutative lumps are not stable dynamically p3[ which leads also to the 
instability of the motion around the moduli space. 

Our approach is free of these drawbacks since we do not make any assumptions 
about the stability. As the analysis shows the dynamics of the system does not look 
as a stable one, moreover, it appears to be stochastic! The regular motion occurs 
only when the distance between lumps is exactly \J 9 In 2. It is interesting to note 
that for some natural initial conditions the dynamics of noncommutative lumps is 
described by finite dimensional matrix model. 

The plan of the actual paper is as follows. First, we introduce the reader to 
the noncommutative lumps in Yang-Mills-Higgs model. After that we analyse the 
lump dynamics in both naive approach when we treat lumps as rigid particles and 
neglect the dynamics of the "shapes" of the lumps and in an exact approach when all 
possible deformations are taken into account. The comparison reveals unexpected 
features in the behaviour of the interacting lumps. We also give directions to 
generalise the description to the case of lumps with arbitrary polarisations. 



2. The Model 

In this paper we consider the noncommutative gauge model described by the 
following action, 

(1) S = j dttrl^X'X' + ^lX'.X^Y 

Here fields X'\ i — 1, . . . , D are time dependent Hermitian operators, acting on 
Hilbert space Ti. which realises a irreducible representation for the onc-dimensional 
Heisenberg algebra generated by, 

(2) [x\x^]^ie. 

Operators x'' satisfying the algebra (^) are said also to be the coordinates of a 
noncommutative two-dimensional plane. In this interpretation the operators of the 
Heisenberg algebra H can be represented through ordinary functions given by their 
Weyl symbols. The composition rule for the symbols is given by the Moyal or star 
product, 

(3) /*5(:r)= e'«'"'^-^^/(x)5(x') 

where /(x) and g{x) are Weyl symbols of some operators, / * g{x) is the Weyl 
symbol of their product and 9^, d'^ denotes the derivatives with respect to x^^ and 
x''^. Integration of a Weyl symbol corresponds to 27r0x trace of the respective 
operator, while the partial derivative derivative with respect to corresponds to 
the commutator, 



(4) 
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where is given by — {l/9)e^^x'^ . Since there is one-to-one correspondence 
between operators and their Weyl symbols we will not distinguish between them, 
i.e. keep the same character for both, unless in the danger of confusion. 

The model (|l|) corresponds to the Hilbert space {N oo limit) of the BFSS 
Matrix Model as well as in different perturbative limits it describes the noncom- 
mutative Yang-Mills(-Higgs) model in the temporal gauge v4o = ^ pO[ . 

Indeed, for equations of motion corresponding to the action ([^), 

(5) X, + ^[X,,[X„X,]]^0, 



one may find static classical solution Xi = Pi, ||2J, |19|, 20 , satisfying, 
(6) [p.„p,]^i9-^, 

with constant invertible 9~j^ . We assume about the solution also the irreducibility 
condition: from [pj, _F] = with all p^, z = 1, . . . , D it follows that _F is a c-number, 
F - I. 

Expanding fields around this solution, Xi = pi + Ai, and Weyl ordering op- 
erators Ai with respect to = 9'^^pj one gets precisely the {D + l)-dimensional 
noncommutative Yang-Mills model for the field given by the Weyl symbol Ai{x). 

Getting another solution with a smaller number of independent pi's, Xa — Pa, 

and Xi = const, i =p+l,p + 2, ~ one gets as a result the model of 

p-dimensional Yang-Mills field interacting with {D — p) scalars. 

Having in mind this equivalence, in what follows we will consider two-dimensional 
form of this noncommutative model. If one forgets for a while also the issues with 
the Gauss law the theory looks like a "simple" noncommutative scalar model in 
(2-1-1) dimensions. 

For our purposes it will be convenient to use two-dimensional "complex coordi- 
nates" given by oscillator rising and lowering operators^ a and a, 

(8) a = -^(xi -l-ix^), a= ^={x^ -ix^), [a,a] = l, 



and the oscillator basis, 



(9) aa \n) — n\n) , a \n) ~ \fn |n — 1) , a |n) = Vn+1 \n + 1) . 

As one can see the solution or (^has divergent traces. Another type of static 
solutions one can find in the model (nl) is given by a configuration with localised 
i.e. lump- like Weyl symbols (in some background Pi)^ It is given by commutative 
matrices of finite ranks [ p^ . Although these lumps carry no energy, — they are 
geometrically nontrivial vacua, we will call them noncommutative lumps due to 
their close relation to ones discussed in the literature [p|-[|l5|. 

Up to a gauge transformation the iV-lump solution is given by, 

N 

(10) X^ = Y.^\\n){n\, 

n=0 



3.4 

of a and a to distinguish them from 



^One should take care that the Gauss law constraints ar e sa tisfied as well. We postpone the 
discussion on the Gauss law constraints until the subsection 

^For the Weyl symbols we will use later z and z instead 
the Hilbert space operators. 

^Fairly speaking these solutions are localised if the fields are treated as scalar ones. Since the 
gauge field definition Aa = Xa — pa implies substraction of a linear function pa this type of 
solutions correspond to functions with linear growth. 
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where cjj is n-th eigenvalue of the (finite rank) operator X*. Due to the finiteness 
of the rank the Weyl symbol of X* vanishes at infinity as quick as Gaussian factor 
times a polynomial. The simplest one-lump solutions can be written in the form, 

(11) Xf'=c,|0)(0|, 

where Ci give the "height" and the "orientation" of the lump, by a proper Lorentz 
transformation, Xi —^ A^X^- it can be made of the form Cj = cSn, while its "polar- 
isation" corresponds to the oscillator vacuum state |0). 

In the star-product form operator ( JIl] ) is represented by the Weyl symbol 

X,{z,z) = 2c,e-2|z|^^ 

The lump shifted along noncommutative plane by a (c-number) vector u is given 
by 

(12) Xf"' = CiC-'f^"" |0) (0| e'f""" = Cie-I"l'e"" |0) (0| e""". 

Its Weyl symbol, correspondingly, is given by X^''\z) = 2c,e-l^-"l'. The shifted 
lump with constant m is a solution again. When u becomes time-dependent one 
can perform the time-dependent gauge transformation to get ^ 

(13) X, ^ e'P^"''Xie-'f''"^ 

which shifts the lump back to the centre, but produce a kinetic term for ^ u'^/2. 
Thus a single noncommutative lump moves freely like a non-relativistic particle. It 
is also stable since its energy at rest is zero. 

In what follows we are going to analyse the situation when there is a couple of 
lumps separated by a distance u. 

3. A PAIR OF INTERACTING LUMPS 

As we have shown in the previous section, a single noncommutative lump can be 
always rotated to have the polarisation |0) and orientation along Xi. When there 
are just two lumps one can choose without loss of generality the configuration to 
involve nontrivially only two matrices e.g. Xi and X2. 

Consider two lumps which are obtained from c|0) (0| by shifts along the non- 
commutative plane by respectively and it(2)- Since the dynamics of the centre 
is free and can be decoupled by a time-dependent gauge transformation similar to 
(|l|), where u = -I- ii(2) is the coordinate of the centre. 

Thus, the configuration we consider looks like, 



(14a) Xi^ciVPV-^ = c\-u/2)(-u/2\, 

(14b) X2^C2V^^PV = c\u/2){u/2\ 

(14c) Xi — const, i — 3, . . . , D, 

where we introduced the shorthand notations, 

(X5) V = e^'/^^P^""^ = g5(afi-au)^ 

(16) P=\0){0\. 



The quotient c can be absorbed by the rescaling of the coupling and the time, 
therefore we can put it generically to unity c = 1 . 

In what follows also Xi i = . . . , D will enter trivially in the equations, so in 
the remaining part of the paper for the simplicity of notations the index i will run 
the range i = 1,2. If we were considering more than two noncommutative lumps 
we would have to keep more matrices. 



^Remember the gauge: Aq = 0. 
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Figure 1. The profile of tlie lump- lump interaction potential in 
the naive approach. 

3.1. Naive picture: rigid lumps. Consider first a naive approach where we 
are dealing with rigid interacting lumps which means that we are neglecting the 
deformations in their shapes. In this case the only parameter which is dynamical 
is the separation distance u. Although, this approximation sounds reasonable, 
later we will consider the exact description which shows that this approach is not 
justified. However, we decided to keep this naive analysis for illustrative purposes. 

To obtain the action describing the dynamics let us insert the ansatz (^4|) into 
the classical action (]l|). The computation of derivatives and traces gives for the 
following M-dependence of the action, 

(17) S[u]^ y'dt(^^H2-2c-^(l-e-^)' 

where we restored the explicit 9 dependence. 

The potential is depicted on the fig.0. According to it sufficiently close lumps 
attract while distant ones repel. At the critical distance Uc = \/20 In 2 they will 
stay in unstable equilibrium. 

The above conclusions concerning the lump dynamics would be valid, however, 
only in the case when one can neglect the involvement of the lump shape in the 
dynamics. To evaluate the importance of the shape dynamics one should consider 
arbitrary deformations of the shape of lumps and separate them from the motion 
of the lump as a whole. 

In the next subsection we analyse the dynamics from the point of view of exact 
field equations of motion. The lump configuration is taken to be the initial condition 
for the field equations. The result we obtain in the next section will invalidate the 
results of the actual naive approach, however, the critical distance Uc will correspond 
to a special case. 

3.2. Exact description: Lumps at rest. The exact description of the lump 
dynamics is given by the field equations of motion for Xi , 

(18) X, + ^[X„[X„X,]]=0 

corresponding to the action (|l|) , supplied with initial conditions given by the lump 
background (p^. Since the equations are second order, in addition to this one has 
to consider the initial values for the time derivatives of Xi. The simplest choice is 
when the configuration at i = is static. Thus, the initial conditions we impose 
are as follows, 

(19a) Xx\t=o^\-u/2){^u/2\, Xi\t=f, = 

(19b) X2\t=o - \u/2) {ul2\ , iTalt^o - 0. 

Considering the lumps in the initial moment as being at rest, produces a con- 
siderable simplification to the equations of motion. Indeed, the initial data dl4) 
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imply that the operators Xt are nonzero only on the two-dimensional subspace 7i„ 
of the Hilbert space which is the linear span of vectors |u/2) and \—u/2). Since, 
in virtue of equations of motion (^8|) , the second time derivative is proportional to 
commutators of Xi then it also vanishes outside the two-dimensional subspace Tiu- 
Due to zero initial conditions for the first derivatives operators Xi will remain all 
the time in the same two-dimensional subspace of the Hilbert space. 

Let us consider only those components of Xi which are nonzero. This reduces 
the Hilbert space operators to ones acting on the two-dimensional subspace Tiu of 
the Hilbert space spanned by |±u/2). Let us introduce an orthonormal basis on 
'Hu- 

The natural orthonormal basis one can build up out of |±it/2), is given by vectors 
|±), defined as follows (see the Appendix), 



(20a) 1+) ^ {]] = ^ {\u/2) + \-ul2)) , 

/2(l + e-5l"l') 

(20b) h) ^ f J) = , ^ i\u/2) - \-u/2)) . 

^2(l-e-5l-l^) 

The singularity in the |— ) in the limit m — s- appears since in this limit \u/2) and 
\—u/2) tends to be parallel and the subspace become one-dimensional. 

In this basis our problem is reformulated in terms of the 2x2 matrix model with 
equations of motion superficially looking the same as 

(21) Xf) + l[xf ,[xf ,xf)]]=0, 

(2) 

but now X^ are finite dimensional 2x2 matrices. The initial conditions in the 
basis (|o|) are rewritten as follows, 

(2) 1 / 1 + e^sl"!^ ^ ^ . . , 

(22a) X^ \t^a = - _ e-l"P ' —liu\^ )' ^i|t=o 0; 




(It is worthwhile to note that the description in terms of 2x2 matrices is valid 
only for the situation where the lumps were initially at rest. Beyond this the 
conditions Xi\t=o = is an additional specification and it says that also the shapes 
of the lumps are not changing at the initial moment. One may try, however relax 
the last condition and consider more general initial configurations including lumps 
with rising/decreasing height for Xi ^ 0. One can consider more general initial 
condition Xilt^Q cx Xi for which the same description in terms of 2x2 matrices 
remains valid. Solutions of this type for u — Q were considered in p5[.) 



1 / l + e-^l"!' Vl -c-l"l' 



Similar equations as ones given by (21) although in a different context of the 



one-dimensional ordinary Yang-Mills model were under study for a long time and 
were initiated by pa-(28|. In the modern context of the application to the finite 



TV matrix model they appear in [E9|-|32 . The system described by such equations 



was shown to exhibit a stochastic behaviour. Let us describe it in more details to 
the application to the present case. 

In order to rewrite the equations (|2^) in the scalar form let us expand the 
matrices Xi in terms of the two-dimensional Pauli matrices CTq,, a = 1, 2, 3, and the 
two-dimensional unity matrix I2 (which in fact is the projector to ?i„) satisfying, 

(23) [o-ajO-p] = le^ii^a^, [fy^h] = 0, 
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The expansion is as follows, 

(24) Xi^2=^lV2+^lV"- 

In terms of this expansion the equations of motion look as follows, 





(25a) 
(25b) 

(25c) 



1 



a 



2 (^"^P 



X2X2f})X^ = 



:{Xl8f,-X1X,p)xfi = 0, 



where XI2 = -'^"2^" 2- For the initial conditions one also has, 
(25d) XZ2 = 0, 



(25e) 

(25f) 
(25g) 
(25h) 



-^^1 |t=0 — -'^2|t=0 



Xl\,=o^~lVl~e-\-\\ 
-^^114=0 — X2\t=o — 0, 



In particular, the equation (25a) says that the scalar parts of the matrices Xi^2 
remains constant during the motion (-'^i 2(^) = 1/2) provided the zero initial con- 
ditions for the "velocities" Xf ~ 0. At the same time the remaining parts are 
subjects to more complicate nonlinear dynamics. 

Before analysing the solutions for 2, a = 1, 2, 3 let us consider their interpreta- 
tion in terms of the lump dynamics over noncommutative space in the star-product 
representation. 

The noncommutative function which corresponds to a particular solution X" (t) 
will be given by 



X,{t;z,z) - il(z,z)+Xf(i)a„(z,z), 



(26) 

where X^lt) are the solutions of to ( |25| ) and I(z, z) with (Tq(x) are the Weyl symbols 
corresponding to the Pauli matrices. 

The respective Weyl symbols are computed in the Appendix. They are given by. 



(27a) 
(27b) 
(27c) 

(27d) 



cri{z,z) 
<^2{z,z) 
(J3{z,z) 



VT^e-l 
2ie-22^ 

Vl-e-l-l^ 



-2N-f|' 



-2|z+f |- 



_ e-2|^- 2 1 + e~2| 

1 - e~l"l V 



-2zz 



+ 



-zu~\-uz\ 



1 



-|u|2 



l{z, z) = aoiz, z) 



1-e- 

2 1 



-2\z+^Y 



2e-2|z|^-i|«| 
1-e-H' 



u—uz zu-\-uz\ 



As it can be seen from the eqs. (^6|), ( |27| ) irrelevant to the particular form of the 
solution Ar"(t), at any time fields Xi{z,z) are nonzero only in the small vicinities 
(of the size of the order of ~ Vd) of of points z = and z = ±u/2. This means 



8 



CORNELIU SOCHICHIU 



that irrespective to the initial distance between the lumps once left with zero initial 
velocities they will not try to leave their places, the dynamics instead will concern 
only their heights and creation of a "baby-lump" in the middle point between them. 
Let us note that this should be surprising as it is in total disagreement with the 
naive approach drawn in the previous subsection, since there is no regime when the 
lumps would behave like rigid particles. 

Let us consider now the the time-dependent functions X"{t) in more details. 
The equations (^) are too complicate to find the general solution, however, for 
our particular initial data one can use the rich symmetry of the model and find a 
simplifying ansatz. 

Assuming that the magnitudes of Xf and X2 are equal Xi{t) = also for 

nonzero times, we can check this assumption later as a consistency condition for 
the ansatz, but also prove it independently of the ansatz using conservation laws, 
one can split X" and into two orthogonal components X" and Y" as follows, 

(28a) = X" + F", X^ = X" - r", 

(28b) X^^{X^+X^), Y=^{X^-X^\ 

the equality of the square modules Xf{t) ~ X2{t) implies that X and Y remain 
orthogonal. The equations of motion in terms of X and Y read, 

(29a) X" = -^Y^X", 

r 

(29b) f"=-4-X^i^", 

r 

where X^ = X°'X" and Y^ = Y°'Y°'. The initial conditions are respectively, 
(29c) X-U=o = ^iXnO) + X°(0)), Y%^o = liXnO) - X^O)), 

(29d) i:"(o) = y"(o) = 0. 

From eqs. ( p9| ) one can see that the directions of X"' and Y°' do not change. 
The fact that X" and are always mutually orthogonal makes the assumption 
X^(t) = -^l(t) for the ansatz ( p8| ) consistent. 

Splitting the vectors X" and Y" in the dynamical magnitude and the static 
direction, given by the constant unimodular vectors 

(30a) =X7^/X^|t^o = (0,0,1), 

(30b) e?.-r"/Vr^|t=o = (1,0,0), 

one has the equations for the magnitudes X and Y, 

(31a) X = -^Y^X, 

9 

(31b) Y = -^X^Y, 

9^ 

which are supplied by the initial data, 



(31c) X|t=o = e-5l-l', Y\t=o = -Vl-e-M\ X\,^o = Yt=o = 0- 

As we mentioned earlier, the system ( |3l| ) exhibits a stochastic behaviour which 
has been studied both numerically and analytically [^-Q. The system is equiv- 
alent to one of a two-dimensional particle in the potential U{X,Y) = X'^Y'^. The 
allowed by energy conservation region of the configuration space is divided into so 
called stadium X ^ Y < 1 where the motion is almost free and four channels along 
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the axes. The motion in channels can be described by the asymptotic formula, in 
the limit when one coordinate is much smaller than another, say X -^Y p3|. 



(32a) 
(32b) 
where 
(32c) 

and 

(32d) 
(32e) 



Y{t) = -^t^ + Wot + Yo 



X{t) 



2A 



9" V ^(0 



+ -^t^Yot + ipo 



A 



^Yi + g^jqYi 

Yo 



^0 




2 

9 ^0^0 



Yo = Yito), 
Wo = Yito), 



to being the time of entrance into the channel. 

In the channel the particle reach the maximal value of y ~ Wq /A after which 
it is reflected back to the stadium. The instability arises when the particle passes 
through the stadium and enters a new channel. So, generally the motion of the 
particle is stochastic. Also there is a discrete set of trajectories which are closed. 
Thus, depending on initial conditions the system can move in a regular periodic 
way, although this motion is unstable as arbitrary small perturbation can push the 
system to the stochastic regime. 

The asymptotic formulae (32a), ( |32b| ) can provide a reliable description of the 
system for a certain period of time for extremal cases when the lump centre sep- 
aration distance is either large {q = e"'"' ^ 1) or small {y/l — <C 1). Thus, 
if u — > oo (q <C 1) then for times less than istoch = gq~^ , one has the asymptotic 
solution. 



(33a) 



(33b) 



Y{t) = 



X{t) 



g2(l-g2)l/2 



X cos 



(q2(l - ^2)1/2/4^2)^2 _ (1 _ ^2)1/2) 



4ff2 



t^'-VT^t 



In the opposite case when the lumps are close one can again give a reliable descrip- 
tion for of the dynamics by the asymptotic formula. 



(34a) 
(34b) 




(q(l - g2)/V)t2 + q 



cos 



9(1-9'), 3 



4g2 



t-" + qt 



valid for times up to of order fgtoch = 5(1 ~ 9^)^^^' after which the system ap- 
proaches the stadium where we cannot control it. 

There is also one particular separation distance which corresponds to periodic 
motion. This happens for the initial conditions X|t=o = —Y\t=o = 1/V2 or u = 
V6'ln2. In this case the motion is periodic and is given by X{t) = Y{t) = f(t) 
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where for f(t) we have the (imphcit) formula, 
(35) fit): • 



Now, let us recall that in terms of X(t) and Y(t) the dynamical field Xi(i, z, z) 
describing the lumps takes according to eq. (|2^) the following form, 

(36a) Xi(i, z, z) = iao(z, z) + X(t)f73(^, z) + r(t)ai(z, z) = 



1-e-^N X(^)- Vl-e-M-y(^) ^^i^^.,. 
l-e-l"l' 
l-e-iN^X(t) + Vl-e-H-y(^) 
^ l-e"l"l' 



1 - e-l"l 
and, 

(36b) X2(i,Z,z) = i(7o(z,z)+X(t)fT3(z,z)-y(t)fTl(z,z) 



1-e-^l-l X(t) + v/l-e-H-y(t) 
l-e-l«l' 



1-e-^N X(0-v/l-e-N-y(t)^_„^_^.,. 



l-e-l"l' 

1 



l_e-l«P 



where the functions 1T1.3 are given by the eqs. (|27|). Let us note that the function 
cri(z, z) is localised in the points where the lumps are i.e. at z = ±m/2 while the 
function (73 (z, z) is localised at both lump positions as well as at the origin where 
is the middle of the lump centres connecting line. 

The analysis of the solution ( |36| ) reveals that once left at their positions the 
lumps will not tend to move away from them but engage in a stochastic change 
of their heights as well as creation of a "baby" lump in the middle point between 
them which is the origin of the noncommutative plane. This process can be reliably 
described for a while of time in the limits when the lumps are placed very close or 
very far, each case degenerating to stochastic, although correlated variation in the 
heights of the lumps. 

3.3. Exact description: Lumps in motion. The difference arising when consid- 
ering moving lumps consists in the initial values for the velocities. Since a generic 
initial condition for the velocities can complicate the system making it back infinite 
dimensional we confine ourselves to such initial configurations which correspond to 
rigid motion of the lumps. 

Thus, one has to replace the initial values for the velocities by the following, 

(37) Xi\t=Q^^^u\t=Q^^^u\t=Q, 



or, explicitly, using (|14| ), 



(38a) Xx\t=ii = -^(■^" + uv)X-Y\t=Q + \^{vaXi_ + vXxa)\t=Q, 

(38b) X2\t=Q = -^(uM + uv)X2\t=Q - \^{vaX2 + wX2a)|t=o, 

where v — u{t = 0), and solve the infinite dimensional operator equation (^). 
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Applying the same strategy as in the case of lumps at rest we see that the initial 
data are given by operators which are nonzero only in a four-dimensional subspace 
of the infinite-dimensional Hilbert space TL, which is spanned by our two old 
vectors |±u/2) together with other two newcomers a |±u/2) = 2{d/du |±u/2)). Let 
us note that they are all linear independent for m 7^ 0. Therefore, the system is 
reduced to the four-dimensional matrix model. 

The treatment of this model differs from what we had with lumps at rest only 
in technical details, therefore we will not discuss it more. 

The qualitative picture one has in this situation does not change much in com- 
parison to the case of lumps a rest. Just as in the previous case there is a stochastic 
dynamics of the heights of the lumps and creation of "baby" -lumps while the centres 
of the lumps will keep moving with the constant velocities. Indeed, for accelerating 
lump operators Xi are nonzero out of the subspace 7i^, which, as we know, does 
not happen. 

In general, the solution is given by a linear combination with time dependent 
coefficients of functions (^l]), their first derivatives {daa/du), (daa/du) and some 
of their second derivatives like {d^aa/ dudu). 



3.4. The Gauss Law. Once we want to relate our system to the the Yang- 
Mills/BFSS model we have to care about the Gauss law constraint, which is ob- 
tained from the variation of the Aq component of the original gauge invariant 
noncommutative Yang-Mills or BFSS action. This constraint looks like follows, 

(39) L=[X,,X,]=Q. 

As we discussed at the beginning of this section the equations of motion imply 
that the quantity ( ^9|) is at least conserved. Indeed, using the equations of motion 
one has 

(40) i = i[X„l,]=0. 

Therefore to get a self-consistent solution for the Yang-Mills/M(atrix) theory 
one has to verify that -/j|t=o = 0- For zero velocity initial conditions this is implied 
automatically, while for the moving lumps one has, 

(41) L = i[X„ Xi\\t^^ = w [(^ ^ a) Xi + + a) X 

where "h.c." stands for the Hermitian conjugate. 

The equation ( ^ ) requires velocities v and v to vanish. However, the violation of 
the Gauss law for nonzero velocities may be interpreted as the presence of nontrivial 
electric charge. Indeed, in the presence of external sources the Gauss law becomes, 

(42) L' = i[X,X,]+p = 0, 

where p is some electric charge density which appears in the action as a term 
A^chargo = / d^^^x pXf^ and which is chosen to cancel ( p] ) exactly. (Here we are 
not going to analyse in which conditions such charge density can be created.) 

As a result we have that the Gauss low is satisfied automatically in the case of the 
lumps at rest, while moving lumps generate some background charge distribution. 



3.5. More dimensions. One can do analogous analysis in more than 2 + 1 dimen- 
sions. 

The only difference which appear in p + 1 dimensions is that one has to compute 
the Weyl symbols of sigma matrices with respect to a different background e.g. one 
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given by (g). As a result one has equations similar to (27) 

2 



(43a) (71 (x) 

(43b) (72 (x) 



(e-(- 



)-G-(x-f ) _ ^-(x+f )-G-(a;+f 



_ 1 

^ 2 

1 _ Q-^u-G-u 



sm u X X, 



9„— M-Gtl , 1 

(43c) (73(x) = ^ f e-(-t)-G-(-t) + e-(^+t)-G-(-+t)) 



COS U X X, 



(43d) 



I(x)=ao(x) = ^— ^^(e 



-(x-f )-G-(x--|) 



+ e" 



-(x+f )-G-(3:+f 



1 _ g-|« G M 



■ COSU X X, 



where we introduced the notation, m x x = 6'^i/u''x'^, ii^v = 1, . . . ,p, and squares 
are computed with the metric G = +V —9^^. In the basis for which the non- 
commutativity matrix 0'^" takes the canonical form, 

\ 

Oo^i(T2 

(44) 6'"' = 

the metric G is diagonal, 

(45) 



Gnu 



^(i)io'2 













6'(2)i(72 













^(3) if 2 












..\ 





6*^25 12 

















V ■■ 









where 
(46) 



1(72 



1 

-1 



and lo 



1 
1 



Then the solution is given by equation similar to ( Pq ) but with I(x) and aa{x) 
instead of respective two-dimensional symbols. 



4. Discussions and Conclusions 

In this paper we considered the dynamics of interacting noncommutative lumps. 

The naive approach for the dynamics is obtained when one considers the motion 
of the lumps as rigid structures and not their deformations. The action in this 
approach is given by initial classical action of the noncommutative model computed 
when all "degrees of freedom" except the positions of the lumps are frozen. In this 
approximation one obtains that the lump pair dynamics is described by a cup- 
shaped potential having minimum at the origin, Gaussian decay at the infinity and 
an unstable equilibrium at the distance y/9hi2. 

The exact analysis of the interacting lump dynamics in the the framework of 
the original noncommutative theory, however, refute above approximation since it 
appears that in fact it is the shape which is affected by the dynamics, while the 
motion of the centres of the lumps is not. 

It is interesting to note that the problem formulated in noncommutative Yang- 
Mills model is reduced to one in finite dimensional matrix model. Thus, in the case 
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of two lumps starting at the rest the exact description reduces to a 2 x 2 matrix 
model. In particular we have that the U(l) part of these model has trivial dynamics, 
while the remaining SU(N) parts generally exhibits stochastic behaviour. 

The property of this dynamics that it does not affect the motion along the line 
connecting the lumps appears to be anti-intuitive to what one could expect from 
interaction of (quasi)particle objects. Let us note that an analogous situation can 
be met in the analysis of the vortex interaction in the applications to solid state 
physics, |36|, where a behaviour similar to one of noncommutative lumps 

was observed long ago. 

It seems that these results can be easily generalised to the case of lumps with 
arbitrary mutual Hilbert space polarisations not related to the shifts along noncom- 
mutative space. The dynamics of such lumps or branes does not differ qualitatively 
from the case shifted ones, however, in this case we do not have simple physical 
picture we can observe. However, from the point of view of mathematical complete- 
ness this would be worthy to be considered, and probably will be done in future 
research. 

It seems that the interpretation in terms of branes when the heights of the lumps 
have the meaning of coordinates of the 0-brane in the direction transversal to the 
noncommutative brane is the most natural. In this context it appears that the dy- 
namics of interacting branes affects only the motion in the transversal directions in 
which DO-branes are stochastically "bouncing" around the noncommutative brane. 

Translating the above said about lumps to the branes, we have learnt that the 
dynamics of two interacting O-branes is described by U(2) M(atrix) model in the 
case when the branes do not move or do not change their polarisations. If the 
branes are affected by motion one needs a matrix model of higher dimension to 
describe it. 

So far we considered only the pair interaction of the lumps. It also would be 
of interest to extend the analysis of the actual paper to a greater number of the 
lumps, eventually to consider the gas of lumps. 
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topoulos, Nikos Tetradis and Ciprian Acatrinei for the friendly atmosphere and 
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discussions on the vortex interactions with Th. Tomaras. I benefited from useful 
discussions at QFT dept of Steklov Mathematical Institute in Moscow. Discussion 
with I.Ya. Aref'eva and P.B. Medvedev improved my understanding the dynamics 
of the SU(2) matrix model. 



Appendix A. Useful formula connecting the two dimensional 

REPRESENTATION WITH OTHER REPRESENTATIONS 

To improve the understanding of the paper we summarise here the formula con- 
necting the three main representations of the objects used in this paper, Hilbert 
space operator, noncommutative functions (Weyl symbols) and two-dimensional 
matrices. 

The two-dimensional space Hu for u 7^ is the span of the two vectors \u/2) = 
e-sl"l Q-^au g^^^ \—u/2) = e~^l"l e^°" |0), where |0) is the oscillator vacuum 
state, to which the lump operators Xi and X2 project to. 

The vectors |±u/2) have unity modules but are not orthogonal. One can eas- 
ily construct an orthonormal basis consisting of vectors {|+) , |— )} given by eqs. 
( |20a| , |20b| ) (see picture |). 
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l-> 



Figure 2. The orthogonal vectors |±) constructed from unimod- 
ular but nonorthogonal |±u/2). 




C7q{z,z) 



0-1 (z,z) 




CT2{Z,Z) 0-3(2, -Z) 

Figure 3. Plots of the profiles of the functions o-o.i,2.3{z, z). 



An arbitrary Hermitian operator acting in this two-dimensional subspace can be 
expanded in terms of ordinary Pauli matrices and unity matrix, 



(47) ^0^1.= (J 'X <^^-r, I 



^2 



-i 

1 



CT3 



1 

-ll ' 



as follows, 

(48) X = X° + X°'aa, 

where X" are computed as, 



(49) 



X" = -trXcr„. 
2 



From the other hand, as operators over the Hilbert space the two-dimensional 
unity matrix^ and Pauli matrices can be expressed as noncommutative functions 
through their Weyl symbols. 



^Which is the projector to the two dimensional subspace of the Hilbert space. 
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The Weyl symbols of operators with bounded square-trace to which undoubtedly 
belong I2 and Ua can be found by a direct formula, 



(50) 



X 



dfcdfce'('^^+''*) trXe-'('^''+'=°^ 



Technically, one can write the matrices in the (nonorthogonal) basis of |±m/2) 
and use the Weyl symbols for the following operators, 



(51) 



|m/2) (u/2| - 2e-2N-"/2|^^ 
\-u/2) (-u/2| - 2e-2|^+"/2|', 
|u/2) (-u/2| ~ 2e-2|^l'+(^"^-^«), 
\-ul2) (u/2| ~ 2e-2|^l'-(^"-^«), 



which can be easily computed. 

The sigma-matrices are expressed in the nonorthogonal basis of |±u/2) as follows, 



(52a) 
(52b) 

(52c) 

and, finally 
(52d) 



1 



0-1 



0^2 



0-3 



-\u[ 



1 



Vl - c-l"l 
e-5l"l 



1 
1 



-.{\u/2) {u/2\ - hV2) {-u/2\), 
^{\u/2) {~u/2\-\^u/2) (V2|), 

(|u/2) (V2| + \-ul2) {-u/2\) + 



c-|«l 

^(|V2)(-V2| + |-V2) (V2|), 



1 



1 - e-l"l 
e-5l-l' 



-(|u/2) {u/2\ + \~u/2) {~u/2\) + 



{\u/2) {^u/2\ + \-u/2) {u/2\) 



1-e-H' 

Inserting ( |5l| ) into (52) one finds immediately the functions (| 
section. 

The plots of functions aa{z, z) can be seen on the fig.§. 
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